BABAR-PUB-05/028 
SLAC-PUB-11332 



Amplitude Analysis of the Decay B ± — > n 7r 

B. Aubert, R. Barate, D. Boutigny, F. Couderc, Y. Karyotakis, J. P. Lccs, V. Poircau, V. Tisserand, and A. Zghiche 
Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, France 

E. Grauges 

IFAE, Universitat Autonoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain 

A. Palano, M. Pappagallo, and A. Pompili 

Universita di Bari, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 

J. C. Chen, N. D. Qi, G. Rong, P. Wang, and Y. S. Zhu 
Institute of High Energy Physics, Beijing 100039, China 

G. Eigen, I. Ofte, and B. Stugu 
University of Bergen, Inst, of Physics, N-5007 Bergen, Norway 

G. S. Abrams, M. Battaglia, A. B. Breon, D. N. Brown, J. Button- Shafer, R. N. Cahn, 
E. Charles, C. T. Day, M. S. Gill, A. V. Gritsan, Y. Groysman, R. G. Jacobsen, R. W. Kadel, 
J. Kadyk, L. T. Kerth, Yu. G. Kolomensky, G. Kukartsev, G. Lynch, L. M. Mir, 

P. J. Oddone, T. J. Orimoto, M. Pripstein, N. A. Roe, M. T. Ronan, and W. A. Wenzel 
Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

M. Barrett, K. E. Ford, T. J. Harrison, A. J. Hart, C. M. Hawkes, S. E. Morgan, and A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

M. Fritsch, K. Goetzen, T. Held, H. Koch, B. Lewandowski, M. Pelizaeus, K. Peters, T. Schroeder, and M. Steinke 
Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44^80 Bochum, Germany 

J. T. Boyd, J. P. Burke, N. Chevalier, and W. N. Cottingham 
University of Bristol, Bristol BS8 1TL, United Kingdom 

T. Cuhadar-Donszehnann, B. G. Fulsom, C. Hearty, N. S. Knecht, T. S. Mattison, and J. A. McKenna 
University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 

A. Khan, P. Kyberd, M. Saleem, and L. Teodorescu 
Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

A. E. Blinov, V. E. Blinov, A. D. Bukin, V. P. Druzhinin, V. B. Golubev, E. A. Kravchenko, 
A. P. Onuchin, S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, and A. N. Yushkov 
Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

D. Best, M. Bondioli, M. Bruinsma, M. Chao, S. Curry, I. Eschrich, D. Kirkby, 
A. J. Lankford, P. Lund, M. Mandelkern, R. K. Mommsen, W. Roethel, and D. P. Stoker 
University of California at Irvine, Irvine, California 92697, USA 

C. Buchanan, B. L. Hartfiel, and A. J. R. Weinstein 
University of California at Los Angeles, Los Angeles, California 90024, USA 

S. D. Foulkes, J. W. Gary, O. Long, B. C. Shen, K. Wang, and L. Zhang 
University of California at Riverside, Riverside, California 92521, USA 

D. del Re, H. K. Hadavand, E. J. Hill, D. B. MacFarlane, H. P. Paar, S. Rahatlou, and V. Sharma 
University of California at San Diego, La Jolla, California 92093, USA 



2 



J. W. Berryhill, C. Campagnari, A. Cunha, B. Dahmcs, T. M. Hong, M. A. Mazur, J. D. Richman, and W. Verkerke 
University of California at Santa Barbara, Santa Barbara, California 93106, USA 

T. W. Beck, A. M. Eisner, C. J. Flacco, C. A. Heusch, J. Kroseberg, W. S. Lockman, G. Nesom, 
T. Schalk, B. A. Schumm, A. Scidcn, P. Spradlin, D. C. Williams, and M. G. Wilson 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

J. Albert, E. Chen, G. P. Dubois-Felsmann, A. Dvoretskii, D. G. Hitlin, 
I. Narsky, T. Piatenko, F. C. Porter, A. Ryd, and A. Samuel 
California Institute of Technology, Pasadena, California 91125, USA 

R. Andrcassen, S. Jayatilleke, G. Mancinelli, B. T. Meadows, and M. D. Sokoloff 

University of Cincinnati, Cincinnati, Ohio 4-5221, USA 

F. Blanc, P. Bloom, S. Chen, W. T. Ford, J. F. Hirschaucr, A. Krciscl, U. Nauenberg, 
A. Olivas, P. Rankin, W. O. Ruddick, J. G. Smith, K. A. Ulmcr, S. R. Wagner, and J. Zhang 
University of Colorado, Boulder, Colorado 80309, USA 

A. Chen, E. A. Eckhart, A. Soffer, W. H. Toki, R. J. Wilson, and Q. Zeng 
Colorado State University, Fort Collins, Colorado 80523, USA 

D. Altenburg, E. Feltresi, A. Hauke, and B. Spaan 
Universitat Dortmund, Institut fur Physik, D-44221 Dortmund, Germany 

T. Brandt, J. Brose, M. Dickopp, V. Klose, H. M. Lacker, R. Nogowski, S. Otto, 
A. Petzold, G. Schott, J. Schubert, K. R. Schubert, R. Schwierz, and J. E. Sundermann 

Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard, G. R. Bonneaud, P. Grenier, S. Schrcnk, Ch. Thiebaux, G. Vasileiadis, and M. Verderi 
Ecole Polytechmque, LLR, F-91128 Palaiseau, France 

D. J. Bard, P. J. Clark, W. Gradl, F. Muheim, S. Playfer, and Y. Xie 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andreotti, V. Azzolini, D. Bettoni, C. Bozzi, R. Calabrese, G. Cibinetto, E. Luppi, M. Negrini, and L. Piemontese 
Universita di Ferrara, Dipartirnento di Fisica and INFN, 1-44100 Ferrara, Italy 

F. Anulli, R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, 
G. Finocchiaro, P. Pattcri, I. M. Pcruzzi,* M. Piccolo, and A. Zallo 
Laboratori Nazionali di Frascati dell'INFN, 1-00044 Frascati, Italy 

A. Buzzo, R. Capra, R. Contri, M. Lo Vetere, M. Macri, M. R. Monge, 

S. Passaggio, C. Patrignani, E. Robutti, A. Santroni, and S. Tosi 
Universita di Genova, Dipartirnento di Fisica and INFN, 1-16146 Genova, Italy 

G. Brandenburg, K. S. Chaisanguanthum, M. Morii, E. Won, and J. Wu 
Harvard University, Cambridge, Massachusetts 02138, USA 

R. S. Dubitzky, U. Langcncgger, J. Marks, S. Schenk, and U. Uwer 
Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

W. Bhimji, D. A. Bowerman, P. D. Dauncey, U. Egede, R. L. Flack, J. R. Gaillard, 
G. W. Morton, J. A. Nash, M. B. Nikolich, G. P. Taylor, and W. P. Vazquez 
Imperial College London, London, SW7 2AZ, United Kingdom 

M. J. Charles, W. F. Mader, U. Mallik, and A. K. Mohapatra 
University of Iowa, Iowa City, Iowa 52242, USA 



3 



J. Cochran, H. B. Crawley, V. Eyges, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin, and J. Yi 

Iowa State University, Ames, Iowa 50011-3160, USA 

N. Arnaud, M. Davier, X. Giroux, G. Grosdidier, A. Hocker, F. Le Diberder, V. Lepcltier, A. M. Lutz, A. Oyanguren, 
T. C. Petersen, M. Pierini, S. Plaszczynski, S. Rodier, P. Roudeau, M. H. Schune, A. Stocchi, and G. Wormser 

Laboratoire de I'Accelerateur Lineaire, F-91898 Orsay, France 

C. H. Cheng, D. J. Lange, M. C. Simani, and D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

A. J. Bevan, C. A. Chavez, I. J. Forster, J. R. Fry, E. Gabathuler, R. Garnet, K. A. George, 
D. E. Hutchcroft, R. J. Parry, D. J. Payne, K. C. Schoficld, and C. Touramanis 
University of Liverpool, Liverpool L69 72E, United Kingdom 

C. M. Cormack, F. Di Lodovico, W. Menges, and R. Sacco 
Queen Mary, University of London, El 4-NS, United Kingdom 

C. L. Brown, G. Cowan, H. U. Flaecher, M. G. Green, D. A. Hopkins, 
P. S. Jackson, T. R. McMahon, S. Ricciardi, and F. Salvatore 
University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

D. Brown and C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 

J. Allison, N. R. Barlow, R. J. Barlow, C. L. Edgar, M. C. Hodgkinson, 
M. P. Kelly, G. D. Lafferty, M. T. Naisbit, and J. C. Williams 

University of Manchester, Manchester M13 9PL, United Kingdom 

C. Chen, W. D. Hulsbcrgcn, A. Jawahery, D. Kovalskyi, C. K. Lac, D. A. Roberts, and G. Simi 
University of Maryland, College Park, Maryland 20742, USA 

G. Blaylock, C. Dallapiccola, S. S. Hertzbach, R. Koflcr, V. B. Koptchev, 
X. Li, T. B. Moore, S. Sarcmi, H. Staengle, and S. Willocq 

University of Massachusetts, Amherst, Massachusetts 01003, USA 

R. Cowan, K. Koeneke, G. Sciolla, S. J. Sekula, M. Spitznagel, F. Taylor, and R. K. Yamamoto 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 

H. Kim, P. M. Patel, and S. H. Robertson 

McGill University, Montreal, Quebec, Canada H3A 2T8 

A. Lazzaro, V. Lombardo, and F. Palombo 

Universita di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 

J. M. Bauer, L. Cremaldi, V. Eschcnburg, R. Godang, R. Kroeger, 
J. Reidy, D. A. Sanders, D. J. Summers, and H. W. Zhao 

University of Mississippi, University, Mississippi 38677, USA 

S. Brunet, D. Cote, P. Taras, and B. Viaud 
Universite de Montreal, Laboratoire Rene J. A. Levesque, Montreal, Quebec, Canada H3C 3,17 

H. Nicholson 

Mount Holyoke College, South Hadley, Massachusetts 01075, USA 

N. Cavallo,''' G. De Nardo, F. Fabozzi,t C. Gatto, L. Lista, D. Monorchio, P. Paolucci, D. Piccolo, and C. Sciacca 
Universita di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 



M. Baak, H. Bultcn, G. Raven, H. L. Snock, and L. Wilden 



4 



NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

C. P. Jessop and J. M. LoSecco 
University of Notre Dame, Notre Dame, Indiana 46556, USA 

T. Allmcndinger, G. Benelli, K. K. Gan, K. Honschcid, D. Hufnagel, P. D. Jackson, 
H. Kagan, R. Kass, T. Pulliam, A. M. Rahimi, R. Ter-Antonyan, and Q. K. Wong 
Ohio State University, Columbus, Ohio 43210, USA 

J. Brau, R. Frcy, O. Igonkina, M. Lu, C. T. Potter, N. B. Sinev, D. Strom, J. Strube, and E. Torrence 

University of Oregon, Eugene, Oregon 97403, USA 

F. Galeazzi, M. Margoni, M. Morandin, M. Posocco, M. Rotondo, F. Simonetto, R. Stroili, and C. Voci 
Universita di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

M. Benayoun, H. Briand, J. Chauvcau, P. David, L. Del Buono, Ch. de la Vaissiere, 
O. Hamon, M. J. J. John, Ph. Leruste, J. Malcles, J. Ocariz, L. Roos, and G. Therin 

Universites Paris VI et VII, Laboratoire de Physique Nucleaire et de Hautes Energies, F-75252 Paris, France 

P. K. Behera, L. Gladney, Q. H. Guo, and J. Panetta 
University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 

M. Biasini, R. Covarelli, S. Pacetti, and M. Pioppi 

Universita di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 

C. Angelini, G. Batignani, S. Bettarini, F. Bucci, G. Calderini, M. Carpinelli, R. Cenci, F. Forti, M. A. Giorgi, 
A. Lusiani, G. Marchiori, M. Morganti, N. Neri, E. Paoloni, M. Rama, G. Rizzo, and J. Walsh 
Universita di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 

M. Hairc, D. Judd, and D. E. Wagoner 
Prairie View A&M University, Prairie View, Texas 77446, USA 

J. Biesiada, N. Danielson, P. Elmer, Y. P. Lau, C. Lu, J. Olsen, A. J. S. Smith, and A. V. Telnov 

Princeton University, Princeton, New Jersey 08544, USA 

F. Bellini, G. Cavoto, A. D'Orazio, E. Di Marco, R. Faccini, F. Ferrarotto, F. Ferroni, M. Gaspero, 
L. Li Gioi, M. A. Mazzoni, S. Morganti, G. Piredda, F. Polci, F. Safai Tehrani, and C. Voena 
Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

H. Schroder, G. Wagner, and R. Waldi 
Universitat Rostock, D-18051 Rostock, Germany 

T. Adye, N. De Groot, B. Franek, G. P. Gopal, E. O. Olaiya, and F. F. Wilson 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

R. Aleksan, S. Emery, A. Gaidot, S. F. Ganzhur, P.-F. Giraud, G. Graziani, G. Hamel de Monchenault, 
W. Kozanecki, M. Legendre, G. W. London, B. Mayer, G. Vasseur, Ch. Yeche, and M. Zito 
DSM/Dapma, CEA/Saclay, F-91191 Gif-sur-Yvette, France 

M. V. Purohit, A. W. Weidemann, J. R. Wilson, and F. X. Yumiceva 
University of South Carolina, Columbia, South Carolina 29208, USA 

T. Abe, M. T. Allen, D. Aston, N. van Bakel, R. Bartoldus, N. Berger, A. M. Boyarski, O. L. Buchmueller, R. Claus, 

J. P. Coleman, M. R. Convery, M. Cristinziani, J. C. Dingfelder, D. Dong, J. Dorfan, D. Dujmic, W. Dunwoodie, 
S. Fan, R. C. Field, T. Glanzman, S. J. Gowdy, T. Hadig, V. Halyo, C. Hast, T. Hryn'ova, W. R. Inncs, M. H. Kelsey, 
P. Kim, M. L. Kocian, D. W. G. S. Leith, J. Libby, S. Luitz, V. Luth, H. L. Lynch, H. Marsiske, R. Messner, 
D. R. Muller, C. P. O'Grady, V. E. Ozcan, A. Perazzo, M. Perl, B. N. Ratcliff, A. Roodman, A. A. Salnikov, 



■5 



R. H. Schindler, J. Schwiening, A. Snyder, J. Stelzer, D. Su, M. K. Sullivan, K. Suzuki, S. Swain, J. M. Thompson, 
J. Va'vra, M. Weaver, W. J. Wisniewski, M. Wittgen, D. H. Wright, A. K. Yarritu, K. Yi, and C. C. Young 

Stanford Linear Accelerator Center, Stanford, California 94309, USA 

P. R. Burchat, A. J. Edwards, S. A. Majewski, B. A. Petersen, and C. Roat 
Stanford University, Stanford, California 94305-4060, USA 

M. Ahmed, S. Ahmed, M. S. Alam, J. A. Ernst, M. A. Saeed, F. R. Wappler, and S. B. Zain 

State University of New York, Albany, New York 12222, USA 

W. Bugg, M. Krishnamurthy, and S. M. Spanier 

University of Tennessee, Knoxville, Tennessee 37996, USA 

R. Eckmann, J. L. Ritchie, A. Satpathy, and R. F. Schwitters 

University of Texas at Austin, Austin, Texas 78712, USA 

J. M. Izen, I. Kitayama, X. C. Lou, and S. Ye 

University of Texas at Dallas, Richardson, Texas 75083, USA 

F. Bianchi, M. Bona, F. Gallo, and D. Gamba 
Universita di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

M. Bomben, L. Bosisio, C. Cartaro, F. Cossutti, G. Delia Ricca, S. Dittongo, S. Grancagnolo, L. Lanceri, and L. Vitale 

Universita di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 

F. Martinez- Vidal 
IFIC, Universitat de Valencia- C SIC, E-46071 Valencia, Spain 

R. S. Panvini* 

Vanderbilt University, Nashville, Tennessee 37235, USA 

Sw. Banerjee, B. Bhuyan, C. M. Brown, D. Fortin, K. Hamano, R. Kowalewski, J. M. Roney, and R. J. Sobie 
University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

J. J. Back, P. F. Harrison, T. E. Latham, and G. B. Mohanty 
Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

H. R. Band, X. Chen, B. Cheng, S. Dasu, M. Datta, A. M. Eichenbaum, K. T. Flood, 
M. Graham, J. J. Hollar, J. R. Johnson, P. E. Kuttcr, H. Li, R. Liu, B. Mellado, A. Mihalyi, 
Y. Pan, R. Prcpost, P. Tan, J. H. von Wimmersperg-Toeller, S. L. Wu, and Z. Yu 

University of Wisconsin, Madison, Wisconsin 53706, USA 

H. Neal 

Yale University, New Haven, Connecticut 06511, USA 
(Dated: August 2, 2005) 

We present a Dalitz-plot analysis of charmless B^ decays to the final state n ± TY ± TY^ using 210 fb _1 
of data recorded by the BABAR experiment at ^/s = 10.58 GeV. We measure the branching fractions 
B{B ± -> tt^tt^) = (16.2 ± 1.2 ± 0.9) x 10~ 6 and B{B ± -> p°(770)-K ± ) = (8.8 ± 1.0 ± 0.6l ( ( ^) x 
10~ 6 . Measurements of branching fractions for the quasi-two-body decays B — > p°(1450)7r ± , 
B ± — > /o(980)7T ± and B ± — > /2(1270)7r ± are also presented. We observe no charge asymmetries 
for the above modes, and there is no evidence for the decays B ± — > XcOTr*, B ± — > /o(1370)7r ± and 

B ± -> CTTT*. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
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I. INTRODUCTION 

The decay of B mesons to a three-body charmless fi- 
nal state offers the possibility of investigating the proper- 
ties of the weak interaction and provides information on 
the complex quark couplings described in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements [1], as well 
as on models of hadronic decays. Measurements of di- 
rect CP-violating asymmetries and constraints on the 
magnitudes and the phases of the CKM matrix ele- 
ments can be obtained from individual decay channels 
in B ± — > 7r ± 7r ± 7r =F [2-5], which are dominated by de- 
cays through intermediate resonances. For example, the 
CKM angle 7 can be extracted from the interference be- 
tween the decay B^ — > XcO 71 "^ which has no CP- violating 
phase, and other modes such as p°(770)tt ± or / (980)7r ± . 
Studies of these decays can also help to clarify the na- 
ture of the resonances involved, not all of which are well 
understood. Of particular interest is whether the a reso- 
nance, which has been observed in other experiments [6- 
8], is also present in B ± — > 7r ± 7r ± 7r T decays. An analysis 
of the full three-body kinematic space is necessary to 
model the interference and extract branching fractions. 

Observations of B- meson decays to the 7r ± 7r ± 7r :F three- 
body final states have already been reported by the Belle 
and BABAR collaborations using a method that treats 
each intermediate decay incoherently [9, 10]. These stud- 
ies have only observed B ± — ► p°(770)n ± , in which other 
possible resonance contributions are treated as back- 
ground. The first measurement of the total branching 
fraction for B^ — > 7r ± 7r ± 7r T was found to be (11 ± 4) x 
10~ 6 [11]. Here, we present results from a full ampli- 
tude analysis for B^ — > 7r ± 7r ± 7r =F decay modes based on 
a 210.3 fb _1 data sample containing 231.6 ± 2.6 million 
BB pairs collected with the BABAR detector [12] at the 
SLAC PEP-II asymmetric-energy e + e~ storage ring [13] 
operating at the Y(4S) resonance at a center-of-mass en- 
ergy of sfs = 10.58 GeV. An additional integrated lu- 
minosity of 21.6 fb _1 was recorded 40 MeV below this 
energy and is used to study backgrounds. The charm 
decay B + — » D a n + , D a — > K + tt~ [14] is used as a cali- 
bration channel as it presents a relatively high branching 
fraction. 



central drift chamber (DCH) in a 1.5-T solenoidal mag- 
netic field. This provides a transverse momentum res- 
olution for the combined tracking system of a PT /pT — 
0.0013pr © 0.0045, where the sum is in quadrature and 
Pt is measured in GeV/ c. 

Charged-particle identification is accomplished by 
combining information on the specific ionization (dE/dx) 
in the two tracking devices and the angle of emission 
of Cherenkov radiation in an internally reflecting ring- 
imaging Cherenkov detector (DIRC) covering the central 
region. The dE/dx resolution from the drift chamber 
is typically about 7.5% for pions. The Cherenkov an- 
gle resolution of the DIRC is measured to be 2.4 mrad, 
for the quartz refractive index of 1.473, which provides 
better than 3 a separation between charged kaons and 
pions over the full kinematic range of this analysis. The 
DIRC is surrounded by an electromagnetic calorimeter 
(EMC), comprising 6580 CsI(Tl) crystals, which is used 
to measure the energies and angular positions of photons 
and electrons. The EMC is used to veto electrons in this 
analysis. 



III. EVENT SELECTION AND 
RECONSTRUCTION 

Hadronic events are selected based on track multiplic- 
ity and event topology. Backgrounds from non-hadronic 
events are reduced by requiring the ratio of Fox- Wolfram 
moments H2/H0 [15] to be less than 0.98. B-meson can- 
didates are reconstructed from events that have four or 
more charged tracks. Each track is required to be well 
measured and originate from the beam spot. They must 
have at least 12 hits in the DCH, a minimum transverse 
momentum of 100 MeV/c, and a distance of closest ap- 
proach to the beam spot of less than 1.5 cm in the trans- 
verse plane and less than 10 cm along the beam axis. 
Charged tracks identified as electrons are rejected. The 
B-meson candidates are formed from three-charged-track 
combinations and particle identification criteria are ap- 
plied. The efficiency of selecting pions is approximately 
95 %, while the probability of misidentifying kaons as pi- 
ons is 15%. The B-meson candidates' energies and mo- 
menta are required to satisfy the kinematic constraints 
detailed in Section V. 



II. THE BABAR DETECTOR 



IV. BACKGROUND SUPPRESSION 



Details of the BABAR detector are described else- 
where [12] . Charged particles are measured with the com- 
bination of a silicon vertex tracker (SVT), which consists 
of five layers of double-sided detectors, and a 40-layer 



tAlso with Universita della Basilicata, Potenza, Italy 
t Deceased 



Backgrounds from e+e~ — ► qq are high and are sup- 
pressed by imposing requirements on event-shape vari- 
ables calculated in the T(4S) rest frame. The first dis- 
criminating variable is cos6*t, the cosine of the angle be- 
tween the thrust axis of the selected B candidate and the 
thrust axis of the rest of the event (all remaining charged 
and neutral candidates). The distribution of |cos6>t| is 
strongly peaked towards unity for qq background whereas 
the distribution is uniform for signal events. We require 
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|cos#t| < 0.65. Additionally, we make requirements on 
a Fisher discriminant T [16] formed using a linear com- 
bination of five variables. The first two variables are 
the momentum-weighted Legendre polynomial moments 
L o = Y,iPi and L 2 = J2iP* x I cos ( i)\ 2 i where p* is the 
momentum of particle i (not from the B candidate) and 
9* is the angle between its momentum and the thrust axis 
of the selected B candidate in the center-of-mass (CM) 
frame. We also use the absolute cosine of the angle be- 
tween the direction of the B and the collision (z) axis in 
the CM frame, as well as the magnitude of the cosine of 
the angle between the B thrust axis and the z axis in the 
CM frame. The last variable is the flavor of the recoiling 
B as reported by a multivariate tagging algorithm [17]. 
The selection requirements placed on |cos#t| and T are 
optimized using Monte Carlo simulated data and have 
a combined signal efficiency of 37% while rejecting over 
98 % of qq background. 

Other backgrounds arise from BB events. The main 
background for our charmless signal events is from charm 
decays, such as three- and four-body B decays involving 
an intermediate D meson, and the charmonium decays 
J/tp -> t + r and ip(2S) -> ( + (r. We remove B candi- 
dates when the invariant mass of the combination of 
any two of its daughters tracks (of opposite charge) is 
within the ranges 3.05 < m n+n - < 3.22GeV/c 2 , 3.68 < 
m, +T - < 3.87GeV/c 2 and 1.70 < m w+n - < 1.93GeV/c 2 , 
which reject the decays J/ip —> C+Cr, tp(2S) -> £ + £" and 
D° — > K + n~ (or 7r + 7r~), respectively. These ranges are 
asymmetric about the nominal masses [18] in order to 
remove decays in which a lepton (I) or kaon has been 
misidentified as a pion. 

We study the remaining backgrounds from charmless 
B decays and from charm decays that escape the ve- 
toes using a large sample of Monte Carlo (MC) simu- 
lated BB decays equivalent to approximately five times 
the integrated luminosity for the data. Any events that 
pass the selection criteria are further studied using ex- 
clusive MC samples to estimate reconstruction efficiency 
and yields. We find that the only significant background 
arises from B — ► K tt^tt decays, in which the kaon 
has been misidentified as a pion. 

We also consider the decay B ± — > K^n^, K® —* tt + it~ 
to be a background, since the K® candidates decay 
weakly and do not interfere with other ir + ir~ resonances 
in B^ — > 7r ± 7r ± 7r =F . We suppress this background by fit- 
ting two oppositely charged pions from each B candidate 
to a common vertex when the invariant mass of the pair 
is below 0.6 GeV/c 2 . This vertex corresponds to the 
decay point for true K® — > tt + -k~ candidates. We remove 
B decays that have fitted Kg candidates with masses be- 
tween 476 and 519 MeV/c 2 (a = 3.6 MeV/c 2 ). 

A further background in this analysis comes from sig- 
nal events that have been misreconstructed by switching 
one or more particles from the decay of the signal B me- 
son with particles from the other B meson in the event. 
The amount of this background is estimated from MC 
studies and is found to be very small; it accounts for 0.6 % 



of the final data sample in the signal region (defined in 
Section V) and is therefore neglected in this analysis. 



V. FINAL DATA SELECTION 

Two kinematic variables are used to select the final 
data sample. The first variable is AE = E* B — y/s/2, 
the difference between the CM energy of the £>-meson 
candidate and y/s/2, where y/s is the total CM en- 
ergy. _Tj3fi_s^cc>nd_J^^^ mass toes 
= y/ (s/2 + Pj • p B ) 2 /Ef - p| where p B is the B mo- 
mentum and (Ei,Pi) is the four-momentum of the ini- 
tial state in the laboratory frame. For signal B de- 
cays, the AE distribution peaks near zero with a res- 
olution of 19 MeV, while the wies distribution peaks 
near the B mass with a resolution of 2.7 MeV/c 2 . The 
mean of the AE distribution is shifted by —5 MeV from 
zero in data as measured from the calibration channel 
B + — > D°n + , D° — ► K + n~ assuming the kaon hypoth- 
esis for the K + candidate. The same shift is also ob- 
served for B+ -f D°ir+, D° -> tt+tt". The typical AE 
separation between modes that differ by substituting a 
kaon for a pion in the final state is 45 MeV, assum- 
ing the pion mass hypothesis. Events in the AE strip 
—65 < AE < 55 MeV are accepted. We also require 
events to lie in the range 5.20 < toes < 5.29 GeV/c 2 . 
This range is used for an extended maximum likelihood 
fit to the toes distribution in order to determine the 
fraction of signal and background events in our data sam- 
ple. The region is further subdivided into two areas: we 
use the sideband region (5.20 < mss < 5.26 GeV/c 2 ) 
to study the background Dalitz-plot distribution and the 
signal region (5.271 < toes < 5.287 GeV/c 2 ) to perform 
the Dalitz-plot analysis. We accept one B-meson can- 
didate per event in the AE strip. Fewer than 3% of 
events have multiple candidates and in those events one 
candidate is randomly accepted to avoid bias. 

The toes signal component is modeled by a Gaus- 
sian function, while the qq background is modeled us- 
ing the ARGUS function [19] with the endpoint fixed 
to the beam energy while the shape parameter is al- 
lowed to float. The BB background toes shape is mod- 
eled with an ARGUS function plus a Gaussian to ac- 
count for the dominant peaking BB background of 86 ±9 
B ± -> K ± n^n ± events, as well as 7 ± 1 B ± -> 7f°7r ± 
events that have K® candidates with invariant masses 
outside the 6 a range. All parameters of the BB compo- 
nent, including the amount of peaking and nonpeaking 
BB background, are obtained and fixed from the MC 
simulation. The fraction of qq events is allowed to float. 
Figure 1 shows the toes projection of the fit to the data 
for B ± — ► 7r ± 7r ± 7r =F . The \ 2 P er degree of freedom for 
this projection is 93/95 and the total number of events 
in the signal region is 1942 (965 and 977 for the B~~ and 
B + samples, respectively). In the signal region, the frac- 
tion of qq background f q q is found to be (71.2 ± 1.8) %, 
while the fraction oiBB backgrounds f BB - is (4.7±0.5) %. 
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FIG. 1: The tties distribution along the AE strip (—65 < 
AE < 55 MeV): the data are the black points with statistical 
error bars, the lower solid red (dark) area is the qq component, 
the middle solid green (light) area is the BB background con- 
tribution, while the upper blue line shows the total fit result. 



The fraction of signal events in the signal region is then 



/si g = l-jW-/ B u = (24.1 ±1.8)% 



VI. DALITZ AMPLITUDE ANALYSIS 

The charmless S-meson decay to the final state 
'k ± 'k ± 'k^ has a number of intermediate states in the 
Dalitz plot [20] that contribute to the total rate, which 
can be represented in the form: 



(IT 



dsi 3 ds 2 3 



\M\ 2 ex 



k 



c k e i6k V k (s 13 ,S2 3 ) 



(1) 



where S13 = to 2 , _ and saa = m 2 , _ are the invariant 
masses squared of the oppositely-charged pion pairs in 
the final state. The invariant mass of each B candidate 
is constrained to the world-average value [18] before S13 
and S23 are calculated. The amplitude for a given decay 
mode k is proportional to c k e l6k 'D k (s 13, S23) with magni- 
tude c k and phase 9 k (— 7r < 9 k < 7r). The distributions 
T> k describe the dynamics of the decay and are a product 
of the invariant mass and angular distributions. For ex- 
ample, if we have a resonance formed from the first and 
third pion from — > 7r ± 7r ± 7r =F , then 



£>fc(si3, S23) = Rk(s 13 ) x T fc (si 3 , S23), 



(2) 



where R k (si3) is the resonance mass distribution and 
T k (si3, S23) is the angular-dependent amplitude. The V k 
are normalized such that 



J |f fc(si3) s 23 )\ 2 dsi 3 ds 2 ?, 



= 1. 



(3) 



The distribution R k (s 13) is taken to be a relativis- 
ts Breit-Wigner lineshape with Blatt-Weisskopf barrier 
factors [21] for all resonances in this analysis except for 
the /o(980), which is modeled with a Flatte lineshape [22] 
to account for its coupled-channel behavior because it 
couples also to the K + K~ channel right at threshold. 
The nonresonant component is assumed to be uniform in 
phase space. The Breit-Wigner function has the form 



R k (s 



1 



13 J 



ml - si 3 - zTO r(s 13 )' 



(4) 



where too is the nominal mass of the resonance and r(si3) 
is the mass-dependent width. In the general case, the 
latter can be modeled as 



r( Sl3 ) = r ( ^ 



2J+1 



rn 



?(<?) 



/si3^(9o) 



(•5) 



The symbol T denotes the nominal width of the reso- 
nance. The values of too and To are obtained from stan- 
dard tables [18]. The value q is the momentum of either 
daughter in the rest frame of the resonance, and is given 
by 

q = - (mj +m§)][si 3 - (ml - m§)]/4si 3 , (6) 

where toi and to 3 are the masses of the two daughter 
particles, respectively. The symbol qo denotes the value 
of q when S13 = m\. The Blatt-Weisskopf barrier pene- 
tration factor Xj(q) depends on the momentum q as well 
as on the spin of the resonance J [21]: 



X a (z) = 1, 

X 1 (z) = V l/(l + z 2 ), 

X 2 (z) = ^l/(z* + 3z 2 



9), 



(7) 

(8) 
(9) 



where z = rq and r is the radius of the barrier, which we 
take to be 4 GeV -1 (equivalent to the approximate size 
of0.8fm). 

In the case of the Flatte lineshape [22] , which is used 
to describe the dynamics of the /o(980) resonance, the 
mass-dependent width is given by the sum of the widths 
in the 7T7r and K K systems: 



r(si 3 ) = r 7r (si 3 ) + T K {s\z), 



(10) 



where 



r7r(si 3 ) = CJttV ' s 13 " 4to 2 , 



Tif(si 3 ) = gx\ s 



''13 



4m 



K 



(11) 



and g v and gx are effective coupling constants, squared, 
for / (980) -> 7T7T and / (980) -> KK, respectively. We 
use the values g n = 0.138 and gn = 4.45(7^ obtained by 
the BES collaboration [8]. 

We use the Zemach tensor formalism [23] for the angu- 
lar distributions of a spin particle (-B ± ) decaying 
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into a spin J resonance and a spin bachelor particle 
(tt ± ). For J = 0, 1, 2, we have [24]: 



1, 



T (0) 

T, (1) = -2p-q, 



.(2) 



(12) 



where p is the momentum of the bachelor particle and 
q is the momentum of the like-sign resonance daughter, 
both measured in the rest frame of the resonance. 

To fit the data in the signal region, we define an un- 
binned likelihood function for each event to have the form 

£(si 3 ,s 23 ) = /sig 

\J2 n k=iCke^V k (s 13 ,s 23 )\\(s 13 ,s 23 ) 
I \J2 k =i c ke i0k T> k (s 13} s 23 )\ 2 e(s 137 s 23 ) ds 13 ds 23 

<2(si3,S23) 



+ /« 

+ Sbb 



'' q J Q(si3,s 23 ) dsi 3 ds 23 
B(s 13 ,s 23 ) 
J B(s 13 , s 23 ) ds 13 ds 23 



(13) 



where n is the total number of resonant and nonrcso- 
nant components in the signal model; e(si 3 ,s 23 ) is the 
signal reconstruction efficiency defined for all points in 
the Dalitz plot; Q(si 3 , s 23 ) is the distribution of qq back- 
ground; B(si 3 , s 23 ) is the distribution of BB background; 
and / s ig, fqq and f B g are the fractions of signal, qq and 
BB backgrounds, respectively. Since we have two iden- 
tical pions in the final state, the dynamical amplitudes, 
signal efficiency and background distributions are sym- 
metrized between s\ 3 and s 23 . The fit is performed al- 
lowing the amplitude magnitudes (c^) and the phases (9i) 
to vary. 

The first term on the right-hand-side in Eq. (13) corre- 
sponds to the signal probability density function (PDF) 
multiplied by the signal fraction / s ; g . This analysis will 
only be sensitive to relative phases and magnitudes, since 
we can always apply a common magnitude scaling factor 
and phase transformation to all terms in the numera- 
tor and denominator of the signal PDF. Therefore, we 
have fixed the magnitude and phase of the most domi- 
nant component, p°(770), to be 1 and 0, respectively. 

As the choice of normalization, phase convention and 
amplitude formalism may not always be the same for 
different experiments, fit fractions are also presented to 
allow a more meaningful comparison of results. The fit 
fraction for resonance k, F kl is defined as the integral of 
a single decay amplitude squared divided by the coherent 
matrix element squared for the complete Dalitz plot as 
shown in Eq. (14). 



F k = 



J\c k e iek V k (s 13 ,s 23 )\ 2 ds 13 ds 23 
1 1 Ej Cje ie Wj(s 13 , s 23 )\ 2 ds 13 ds 23 ' 



(14) 



where the integrals arc performed over the full kinematic 
range. Note that the sum of these fit fractions is not nec- 
essarily unity due to the potential presence of net con- 
structive or destructive interference. 



VII. DALITZ PLOT BACKGROUNDS AND 
EFFICIENCY 

The dominant source of background for this analysis 
comes from qq events. We use a combination of on- 
resonance sideband data and off-resonance data to get 
the background distribution for the Dalitz plot. Note 
that for the on-resonance sideband data, we subtract any 
contributions from BB background (from MC), since this 
is handled separately. Since the background peaks at the 
edges of the Dalitz plot, we use a coordinate transfor- 
mation to a square Dalitz plot in order to improve the 
modeling of the background distribution. Considering 
the decay B + — > 7r + 7r + 7r~, the new coordinates are m' 
and 9' , which arc defined as 



m 



9' 



1 

— cos 

7T 



i m ++ — m ++ [mm| 
m ++ [max] — m ++ [min] 



1 



(15) 



where m ++ is the invariant mass of the like-sign pions, 
m ++ [max] = raj — m^ and m ++ [min] = 2m T are the 
boundaries of m ++ , while 9 ++ is the helicity angle be- 
tween the momentum of one of the like-sign (ir + ) pions 
and the ir~ momentum in the 7T+7T+ rest-frame. Note 
that the new variables range from to 1. The Jacobian 
transformation J between the normal Dalitz plot vari- 
ables to the new coordinates is defined as 



dsi 3 ds 23 = \ J\dm'd9' . 
The determinant \J\ of the Jacobian is given by 



l-^l =4|Pi|jj^|m ++ - 



dm ++ dcos9 ++ 



dm' 



89' 



(16) 



(17) 



where \p\ | is the momentum of one of the 7r + candidates 
and \p 2 \ is the momentum of the tt~ track, both mea- 
sured in the rest frame of the tt + tt + system. The partial 
derivatives in Eq. (17) are given by 



dm 



++ 



dm' 
d cos 9 ++ 

d9' 



— — sin("7r77j')(m ++ [max] — m ++ [min]), 



= — 7rsin("7r6>'). 



(18) 



We get similar expressions for B~ — > 7r~7r~7r + . Fig- 
ure 2 shows the qq background distribution, obtained 
by combining on-resonance sideband and off-resonance 
data. Figure 3 shows the BB background distributions, 
which originate from B^ — > _ft' ± 7r =F 7r ± and B^ — > K^tt^ 
decays. Note that the peaks along the edges of the nor- 
mal Dalitz plot distribution are more spread out in the 
square Dalitz plot format. We use the latter to repre- 
sent the qq and BB backgrounds in the amplitude fit, 
applying linear interpolation between bins. 

The signal efficiency e(s 13 , s 23 ) used in Eq. (13) is mod- 
eled using a two-dimensional histogram with bins of size 
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FIG. 2: Dalitz plot of the qq background obtained from on-resonance sideband and off-resonance data. The left plot shows the 
distribution in normal Dalitz plot coordinates, while the right plot shows the equivalent distribution in the new square Dalitz 
plot coordinates, defined in Eq. (f5). The empty regions correspond to events removed by the charm vetoes. The area of each 
small square is proportional to the number of events in that bin. 
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FIG. 3: Dalitz plot of the BB background obtained from Monte Carlo simulated events. The left plot shows the distribution in 
normal Dalitz plot coordinates, while the right plot shows the equivalent distribution in the new square Dalitz plot coordinates, 
defined in Eq. (15). The empty regions correspond to events removed by the charm vetoes. The area of each small square is 
proportional to the number of events in that bin. 



0.4 ( GcV/c 2 ) 2 x 0.4 ( GeV/c 2 ) 2 and is obtained using 1.1 lection criteria are applied except for those corresponding 
million B ± — > 7r ± 7r ± 7r =F nonresonant MC events. All se- to the invariant-mass veto regions mentioned in Sec. IV. 
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The efficiency at a given bin is defined as the ratio of the 
number of events reconstructed to the number of events 
generated in that bin. Corrections for differences between 
MC and data in the particle identification and tracking 
efficiencies are applied. The efficiency shows little vari- 
ation across the majority of the Dalitz plot, in which 
the average efficiency is measured to be (13.00 ± 0.04) %, 
however there are decreases towards the corners where 
one of the particles has a low momentum. The effect of 
experimental resolution on the signal model is neglected 
since the resonances under consideration are sufficiently 
broad. No difference in efficiency is seen between B~ and 
B + decays at the 2% level. 



VIII. PHYSICS RESULTS 

We fit the B~ and B + samples independently to ex- 
tract the magnitudes and phases of the resonant and non- 
resonant contributions to the charmless B ± — > 7r ± 7r ± 7r T 
Dalitz plot, using Eq. (13). The nominal fit model con- 
tains the resonances p°(770), p°(1450), /o(980), / 2 (1270) 
and a uniform nonresonant contribution. This is chosen 
using information from established resonance states [18] 
and the x 2 variation observed when omitting one of the 
five components. The x 2 value is calculated using the 
formula 

2 v 5 ^ bfi - Z(^)] 2 

1=1 Gi 

where yi is the number of events in bin i of the invariant 
mass or Dalitz-plot distribution, f{x{) is the expected 
number of events in that bin as predicted by the fit re- 
sult and Oi is the error on yi {y/yl). The number of 
degrees of freedom (nDof) is calculated as rif, — k — 1, 
where n& is the total number of bins used and k is the 
number of free parameters in the fit (4 magnitudes and 
4 phases). A minimum of 10 entries in each bin is re- 
quired; if this requirement is not met then neighboring 
bins are combined. Typically, n& is equal to 35 and 75 for 
the invariant mass and Dalitz-plot distributions, respec- 
tively. Since we observe no charge asymmetry in the qq 
and BB backgrounds, we use the charge-averaged back- 
ground distributions shown in Figs. 2 and 3 for the B~ 
and B + fits. The results of the nominal fit to B~ and 
B + on-rcsonance data in the signal region are shown sep- 
arately in Table I. From Eq. (14), it can be seen that the 
fit fraction statistical uncertainty will not only depend 
on the uncertainties of the magnitude and phase of the 
given resonance, but also on the statistical errors of all 
amplitudes. Therefore, we use a MC pseudo-experiment 
technique to obtain the statistical uncertainty on each 
fit fraction. Each pseudo-experiment is a sample of MC 
generated events that contains the correct mixture of sig- 
nal and background, which are distributed across the 
Dalitz plot according to the PDFs defined in Eq. (13). 
We fit these MC samples and plot the distributions of fit 



fractions Fk obtained from a thousand such experiments. 
The statistical uncertainty for each Fk is then the value 
of the width of the Gaussian function that is fitted to the 
Fk distribution. 

Figure 4 shows the mass projection plots for the 
nominal fits to B~ and B + data, while Fig. 5 shows 
the background-subtracted Dalitz plot of the combined 
B^ — ► 7T^~ 7T^~ tt^ 1 data in the signal region. The x 2 /nDof 
values for the opposite-sign and like-sign invariant mass 
projections for B~ (B+) are 51/34 and 27/37 (35/35 and 
47/35), respectively. The x 2 / n Dof values for the two- 
dimensional Dalitz plots are 74/74 and 70/75 for B~ and 
B + , respectively. The four resonant contributions plus 
the single uniform phase-space nonresonant model are 
able to describe the data adequately within the statisti- 
cal uncertainties. For a given resonance, the comparison 
of the fit fraction, not the magnitude, to its uncertainty 
gives a measure of how significant its contribution is to 
the Dalitz plot. Note that the fit fraction uncertainties 
shown in Table I are larger than the uncertainties of the 
magnitudes. This is due to the dependence of the former 
on all of the other amplitudes, via the denominator in 
Eq. (14). It can be clearly seen that the dominant con- 
tribution to the charmless B ± — ► 7r ± 7r ± 7r T Dalitz plot is 
from the p°(770) resonance. Approximately 10% of the 
p°(770) fit fraction lies in the tail of the mass distribution, 
defined as the region outside tuq ± 3Fo. In addition, the 
fraction of p°(770) within one width of the / 2 (1270) res- 
onance line-shape is approximately 13%, which is equiv- 
alent to half of the /2(1270) fit fraction. Further fits are 
performed to the data by removing one two-body com- 
ponent at a time from the nominal model. Removing 
the p°(770), /2(1270) or nonresonant components give 
significantly poorer fit results. Omitting the p°(1450) or 
/o(980) components, which are present at the 1.5 a level, 
gives a small change in the goodness-of-fit x 2 (Eq. (19)). 

We have also tested the introduction of the Xco & n d 
/o(1370) resonances, as well as the low-mass tt + it~ pole, 
known as the a. Analysis of data from the E791 ex- 
periment for D + — > 7r + 7r~7r + [7] and recent data from 
the BES collaboration for J/tp — > cj7r + 7r~ [8] show evi- 
dence of the a. Also a large concentration of events in the 
1 = S-wave irir channel has been seen in the region 
around 500 — 600 MeV in pp collisions [6] . This pole is 
predicted from models based on chiral perturbation the- 
ory [25], in which the resonance parameters are M — 
iT/2 = [(470 ±30) -i(295±20)] MeV. Consequently, 
the a resonance is predicted in B ± — > n^n^TT^ decays. 
For this Dalitz-plot analysis the a resonance is modeled 
using the parameterization suggested by Bugg [26] . The 
contributions that these three resonances make to the 
nominal fit results are not significant and so we place 
upper limits on them. 

To make comparisons with previous measurements and 
theoretical predictions it is necessary to convert the fit 
fractions into branching fractions. These are estimated 
by multiplying each fit fraction by the total branching 
fraction for the B~ and B + fits, which arc then aver- 
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TABLE I: Results of the nominal fits to B — > n n n and B + — > 7r + 7r + 7r data. The first errors are statistical, while the 
second and third uncertainties are systematic and model-dependent, respectively, all of which are detailed in Sec. IX. All phases 
are in radians. 



Component B Fit Result B + Fit Result 



1.0 

7.9 



p (770) Fraction (%) 50.6 ± 7.3 ± 2.2±£* 57.8 ± 6.8 ± 3.5^ 

p°(770) Magnitude 1.0 (fixed) 1.0 (fixed) 

p°(770) Phase 0.0 (fixed) 0.0 (fixed) 

p°(1450) Fraction (%) 6.8 ± 4.5 ± 1.8 ± 0.6 4.9 ± 5.5 ± 1.4±§;| 

p°(1450) Magnitude 0.37 ± 0.11 ± 0.05 ± 0.02 0.29 ± 0.17 ± 0.06 ± 0.08 

p°(1450) Phase +1.99 ± 0.57 ± 0.10 ± 0.08 +0.31 ± 0.70 ± 0.15 ± 0.34 

/o(980) Fraction (%) 3.8 ± 4.9 ± 0.9 ± 2.1 11.1 ± 5.1 ± ;^ 

/o(980) Magnitude 0.27 ± 0.10 ± 0.05 ± 0.07 0.44 ± 0.12 ± 0.03 ± 0.08 

g;^ -0.79 + 0.62 + 0.171^ 



/o(980) Phase -1.59 + 0.47 + nB+0 - 15 



/ 2 (1270) Fraction (%) 14.2 ± 4.6 ± 1.3 ± 0.5 14.1 ± 4.8 ± 1.4±g;X 

/ 2 (1270) Magnitude 0.53 ± 0.10 ± 0.02 ± 0.03 0.49 ± 0.11 ± 0.02 ± 0.05 

/ 2 (1270) Phase +1.39 ± 0.41 ± 0.09 ± 0.09 +1.85 ± 0.47 ± 0.12±g;g^ 

Nonrcsonant Fraction (%) 15.0 ± 8.6 ± 1.9±t;l 12 - 6 ± 7 - 1 ± 2 - G -H 

Nonresonant Magnitude 0.54 ± 0.13 ± 0.03 ± 0.09 0.47 ± 0.14 ± 0.05 ± 0.06 

Nonresonant Phase -0.84 ± 0.38 ± 0.06 ± 0.04 -2.80 ± 0.46 ± 0.07 ± 0.07 



aged. The total branching fractions B^ ot and B^ ot for 
B~ — > 7r~7r~7r + and B + — ► 7r + 7r + 7r~, respectively, are 
defined as 



B. 



Slg 



(20) 



where A 7 "* is the total number of events in the signal 
region, / s i g is the signal fraction defined earlier, N B g is 
half the total number of BB pairs in the data sample [27] 
and (e ± ) is the average efficiency across the Dalitz plot 
weighted by the fitted signal distribution, which is equal 
to (12.4 ± 0.1)% for B~ and B+ . The average total 
branching fraction is then just equal to j{B^ ot + B^ ot ), 
while the average branching fraction for each resonance 
k is given by 



B k = l(F-B^+F+B+ t ), 



(21) 



where F k (F^) is the fit fraction for resonance k for B 
(B+). 

For components that do not have statistically signifi- 
cant fit fractions 90% confidence-level upper limits are 
evaluated. Upper limits are also found for the xm, 
/o(1370) and a components. These limits are calculated 
by generating many pseudo-MC experiments from the re- 
sults of fits to the data, with all systematic sources (see 
Sec. IX) varied within their 1 a Gaussian uncertainties. 
We fit these MC samples and plot the fit fraction distribu- 
tions. The 90% confidence-level upper limit for each fit 
fraction is then that which removes 90 % of the pseudo- 
MC experiments. A branching fraction upper limit is 
then the product of the upper limit on a fit fraction with 



the total branching fraction B^ ot . Corrections applied to 
the signal efficiency due to differences between data and 
MC are described in Sec. IX. We include the variation 
of (e^) due to these corrections by using another large 
set of pseudo-MC experiments, which is generated and 
fitted to the Dalitz-plot model. The content of each bin 
in the efficiency histogram is increased (decreased) by 
the same random fluctuation given by the uncertainty of 
the efficiency correction (5.1%). The 90% confidence- 
level upper limit on the value of the reciprocal of the 
efficiency (1/0.117) is taken as the value of l/(e ± ) for 
the total branching fraction calculation given in Eq. (20) 
that is then used to find the upper limits for the res- 
onance branching fractions. If the upper limits differ 
between B~ and B + , we choose the larger value to be 
conservative. 

In addition to fit fractions and phases, the charge (CP) 
asymmetries for the signal model components are also 
measured. The charge asymmetry for the total branching 
fraction is defined as 



A 



N 



N + 

sig 



sig J "sig 



(22) 



where A^ sig (A^t g ) is the number of signal events for the 
B~ (B + ) sample. The charge asymmetries for the fit 
fractions arc defined as 



A k = 



Frm 



F^N^r 

k sig 



F k ~N^ + F+Nf s 



(23) 



The measured branching fractions and charge asymme- 
tries are summarized in Table II. The total branch 
ing fraction of the charmless i? ± —* 7r ± Tr ± 7r = F 



7r ± 7r ± 7r =F decay, 
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FIG. 4: Projection plots of the fit results for B~ — * h~h~ti and B — > ir it ir~ onto the mass variables m n + n - and m n - w - 
(m^.+^+). The upper (lower) plots are for the B~ (B + ) data sample. The data are the black points with statistical error 
bars, the lower solid red (dark) histogram is the qq component, the middle solid green (light) histogram is the BB background 
contribution, while the upper blue histogram shows the total fit result. The large dips in the spectra correspond to the charm 
vetoes. 



(16.2 ± 1.2 ± 0.9) x 1(T 6 , is consistent with the cur- 
rent world-average value of (11 ± 4) x 10~ 6 [18]. The 
measured branching fraction for the decay B — > 
p°(770)7r ± ,p°(770) -► 7T+7T-, (8.8±1.0±0.6±^)x lO" 6 , 
agrees with the world-average value of (8.6 ± 2.0) x 



10 -6 [18] and is consistent with the average theoretical 
predictions of 11. 9x 10~ 6 and 8.4x 10 -6 that are based on 
QCD factorization [28] and pole-dominance models [29], 
respectively. The upper limits reported for the other res- 
onance modes are an order of magnitude lower than pub- 
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FIG. 5: Background subtracted Dalitz plot of the combined 
B — > 7r ± 7r ± 7r =F data sample in the signal region. The Dalitz 
plot is symmetrized about the y = x axis. The empty regions 
correspond to events removed by the charm vetoes. 



lished limits [18]. The total charge asymmetry has been 
measured to be consistent with zero to a higher degree 
of accuracy than previous measurements [11]. A repre- 
sentative theoretical value of the charge asymmetry for 
B^ — ► p°(770)7r ± is +4.1% [28], ignoring uncertainties 
due to weak annihilation processes, in agreement with 
our measurement. 



IX. SYSTEMATIC STUDIES 

The systematic uncertainties that affect the measured 
fit fractions, amplitude magnitudes and phases are eval- 
uated separately for B~ and B + . The first source of 
systematic uncertainty is the modeling of the signal ef- 
ficiency. The charged-particle tracking and particle- 
identification fractional uncertainties are 2.4 % and 4.2 %, 
respectively. The first is estimated by finding the dif- 
ference between data and MC of the track-finding ef- 
ficiency of the DCH from multihadron events. A pre- 
cise determination of the DCH efficiency can be made 
by observing the fraction of tracks in the SVT that are 
also found in the DCH. The probabilities of identify- 
ing kaons and pions is measured using the decay mode 
D* + —> D°n + , D° —> K~n + , which provides a very pure 
sample of pions and kaons. The difference observed be- 
tween data and MC for the kaon and pion efficiencies 
gives the combined systematic uncertainty of 4.2 % for 
our signal mode. There are also global systematic errors 
in the efficiencies due to the criteria applied to the event- 
shape variables (1.0%) and to AE and tues (1.0%). The 
total fractional systematic uncertainty for the efficiency 
from these sources is 5.1%. Corrections due to differ- 
ences between data and MC have also been included for 
the selection requirements on cos 6t, J 7 , AE and Toes- 



These are found by comparing the difference in the se- 
lection efficiency between data and MC for the control 
sample B+ -> D°tt+. 

The variation of the efficiency across the Dalitz plot is 
also evaluated by performing a series of fits to the data 
where the efficiency histogram has each bin fluctuate in 
accordance with its binomial error. This introduces an 
absolute uncertainty of 0.01 for the magnitudes, 0.02 to 
0.05 for the phases, and a fractional uncertainty between 
1 % and 4 % for the fit fractions. For the average effi- 
ciency, and hence for the total branching fraction, this is 
a very small effect, evaluated at 0.1 %. 

The next source of systematic uncertainty comes from 
the modeling of the backgrounds. The systematic uncer- 
tainty introduced by the BB background and qq back- 
ground has two components, each of which can poten- 
tially affect the fitted magnitudes and phases differently. 
The first component arises from the uncertainty in the 
overall normalization of these backgrounds, while the sec- 
ond component arises from the uncertainty on the shapes 
of the background distributions in the Dalitz plot. The 
uncertainties on the magnitudes, phases and fit fractions 
due to the normalization uncertainty are estimated by 
varying the measured background fractions in the signal 
region by their statistical errors. The maximum uncer- 
tainty for the magnitude (phase) is 0.03 (0.02) due to 
the qq background normalization uncertainty and 0.01 
(0.01) due to the BB background normalization uncer- 
tainty. These uncertainties are added in quadrature. The 
fit fractions have relative uncertainties in the range 1 % to 

9 %. The uncertainties on the fit fractions and phases due 
to the Dalitz-plot background distribution uncertainty is 
estimated in the same way as the efficiency variation, 
namely varying the contents of the histogram bins in ac- 
cordance with their Poisson errors. To be conservative, 
each magnitude (phase) has been given an uncertainty of 
0.02 (0.02) due to the qq background distribution uncer- 
tainty and 0.02 (0.01) due to the BB background distri- 
bution uncertainty, which are then added in quadrature. 
The fit fractions have relative uncertainties ranging from 
1% to 10%. 

To confirm the fitting procedure, 1000 MC pseudo- 
experiments are created from the fitted magnitudes and 
phases and each sample is fitted 100 times with random- 
ized starting parameters. A fit bias of approximately 

10 % is observed for some of the smaller components and 
is included in the systematic uncertainties for the mag- 
nitudes, phases and fit fractions. 

There is a range of different values for the coupling con- 
stants and gx for the Flatte description of the / (980) 
resonance [8, 30, 31]. A model-dependent systematic un- 
certainty is assigned for all magnitudes, phases and fit 
fractions based on the differences between the results of 
the nominal fit and those when the different coupling con- 
stants for the /o(980) are used. There is also the question 
of whether the nonresonant component has an amplitude 
that varies across the Dalitz plot. For the nominal fit, 
uniform phase-space is used for this component in the 
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TABLE II: Summary of average branching fraction (£>) and charge asymmetry (.4) results. The first uncertainty is statistical, 
the second is systematic, while the third is model-dependent. 



Mode 






23(5* -> Mode)(HT 6 ) 


90% CL UL B (10- 6 ) 


A (%) 


B ± -> Tv ± n ± n^ Total 






16.2 ± 1.2 ±0.9 




-0.7 ±7.7 ±2.5 


p°(770)7r ± ,p°(770) 7v- 






q o _i_ i n -t- n ft+o- 1 
o.o m i.u m u.o_q ^ 






p° ( 1450) Tr^p 1 ^ 1450) -> 






1.0 ±0.6 ±0.2 ±0.2 


< 2.3 


+15.5 ± 62.1 ±7.9t?;o 


/o(980) 7 r ± ,/ (980)^7r- 






1.2 ±0.6 ±0.1 ±0.4 


< 3.0 


-49.5 ± 53.7 ±4.9t|;£ 


h (1270)^, h (1270) 






2.3 ±0.6 ±0.2 ±0.3 


< 3.5 


-0.4 ± 24.7 ±2.8l?;g 


S — ^ ^-k^-k^ 1 Nonresonant 


2.3 ±0.9 ±0.3 ±0.4 


< 4.6 


±8.0 ±41.2 ±6.5 ±2.4 


XcO^^cO -» 7T + 7T 








< 0.3 




/o(1370)7T ± ,/ (1370) 




7T 




< 3.0 




<77T , (T — > 7T + 7T~ 








< 4.1 





absence of any a priori model. An alternative parame- 
terization gives the nonresonant dynamical amplitude to 
be of the form 

V NR (s 13 , s 23 ) = e- QS13 + e' as -\ (24) 

where a is a constant [32]. This parameterization does 
not give significant differences compared to the nominal 
fit results (Table I) for a = 0.11 ± 0.02, which is the 
average of the values shown in Ref. [32] . These differences 
are included in the model-dependent systematic error, as 
well as when the \c0, /o(1370) or a resonances are added 
to the fit. 

The dominant systematic uncertainty for the total 
branching fraction B^ ot is due to the efficiency corrections 
(5.1 %). There is a 1 % fractional error on the weighted 
efficiency (e) due to the statistical uncertainties of the 
fitted amplitudes of the various components. There is an 
additional uncertainty in the value of N BB , evaluated at 
1.1 %, as well as the fractional uncertainty in the amount 
of BB background present (2%). The systematic uncer- 
tainties for the resonance branching fractions B^ are just 
the quadratic sum of the systematic errors for the res- 
onance fit fractions Fk and all the contributions to the 
systematic error for the total branching fraction B tot ex- 
cept the fixed BB background component, since this is 
already included in the fit fraction systematics. 

For the charge asymmetries, systematic uncertainties 
from fit biases, efficiency corrections and fluctuations 
in the background and efficiency histograms are not in- 
cluded, since they cancel out. Finally, an uncertainty of 
2 % is assigned for the total and fit fraction charge asym- 
metries due to a possible detector charge bias, which is 
determined by finding the difference between the total 
number of positively and negatively charged tracks in 
the on-resonance data sample. 

X. SUMMARY 

The total branching fraction for the charmless decay 
B ± -> TT ± n ± n^ is measured to be (16.2±1.2±0.9) x 10~ 6 , 



where the first uncertainty is statistical and the second is 
systematic. The dominant component in the charmless 
B ± — > 7T^~ 7T^~ tt^ 1 Dalitz plot is the p°(770) resonance. We 
have a 3 (J indication for the presence of the /2(1270) 
and nonresonant components. The fit fractions of the 
resonances p°(1450) and /o(980) are not statistically sig- 
nificant. The decay B ± — > p°(770)n ± has a measured 
branching fraction of (8.8 ± 1.0 ± 0.6^'j) x 10~ 6 , which 
is consistent with previous measurements [9, 10] and the- 
oretical calculations [28, 29]. This decay can be used to 
help reduce the theoretical uncertainties in the extraction 
of the CKM angle a from the neutral decays B° —* p ±r K^ 
and B° — > p°7r° [5]. It is found that there is no con- 
tribution from the Xco resonance to the — > n ± n ± TT T 
Dalitz plot, which means that the methods advocated in 
Ref. [2-4] to measure the CKM angle 7 are not feasible 
with our current dataset. There is also little evidence 
for contributions from the /o(1370) and a resonances. 
Differences in the parameterizations of the /o(980) and 
nonresonant components do not significantly affect the 
results. Charge asymmetries observed for the total rate 
and resonance fit fractions are consistent with zero, and 
90 % confidence-level upper limits are provided for the 
branching fractions for resonances that do not have sta- 
tistically significant fit fractions. The results presented in 
this paper supersede those of previous BABAR analyses. 
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